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Abstract 

We determine the distribution of cluster sizes that emerges from an initial phase of homogeneous 
aggregation with conserved total particle density. The physical ingredients behind the predictions 
are essentially classical: Super-critical nuclei are created at the Zeldovich rate, and before the 
depletion of monomers is significant, the characteristic cluster size is so large that the clusters 
undergo diffusion limited growth. Mathematically, the distribution of cluster sizes satisfies an 
advection PDE in "size-space". During this creation phase, clusters are nucleated and then grow 
to a size much larger than the critical size, so nucleation of super-critical clusters at the Zeldovich 
rate is represented by an effective boundary condition at zero size. The advection PDE subject 
to the effective boundary condition constitutes a "creation signaling problem" for the evolving 
distribution of cluster sizes during the creation era. 

Dominant balance arguments applied to the advection signaling problem show that the charac- 
teristic time and cluster size of the creation era are exponentially large in the initial free-energy 

2 G* 

barrier against nucleation, G*. Specifically, the characteristic time is proportional to e 5fc s T and 
the characteristic number of monomers in a cluster is proportional to e 5fc s T . The exponentially 
large characteristic time and cluster size give a-posteriori validation of the mathematical signaling 
problem. In a short note, Marchenko [4| obtained these exponentials and the numerical pre-factors, 
| and |. Our work adds the actual solution of the kinetic model implied by these scalings, and 
the basis for connection to subsequent stages of the aggregation process after the creation era. 
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I. INTRODUCTION 



The traditional idea of nucleation is the growth of clusters in an inexhaustible monomer 
bath, by fluctuations over a high free-energy barrier G*. The inexhaustibly of the monomer 
bath means that the super-saturation — and hence G> — are constant in time. In the limit 



G^/ksT 3> 1, Zeldovich 13) derived an asymptotic result for the steady nucleation rate per 



unit volume, proportional to e k s T . The exponential smallness of the nucleation rate in 
is the a-posteriori justification of the analysis. 

The steady Zeldovich rate is only a first step towards a large picture. Starting from pure 
monomer, there is a duration of so-called "transient nucleation" in which the nucleation 
rate ramps up from zero to the steady Zeldovich value. Asymptotic analysis of transient 
nucleation was carried out by Shneidman Q and then by Bonilla et al Q. This is the very 
"beginning" , but what is the "end" of nucleation in a closed system? 

Qualitatively, it is clear: The super-saturation decreases due to the growth of clusters, 
thereby increasing the free-energy barrier, and eventually the nucleation of new clusters is 
shut down. We refer to the duration between the aforementioned transient nucleation and 
this "shutting down" as the "creation era." We use the word "creation" since conventionally 
"nucleation" refers to the steady state process in an inexhaustible bath as first studied by 
Zeldovich. 

This paper proposes a quantitative theory of the creation era. It is quasi-static: We 
assume that the Zeldovich rate applies even though the super-saturation slowly decreases in 
time, and the free energy barrier slowly increases. The quasi-static assumption is affirmed 
a-posteriori using Shneidman's criterion [h]] in appendix [B] Since the nucleation rate is the 
exponential of a large negative quantity — a small relative increase in is sufficient 
to reduce the nucleation rate to a small fraction of its initial value. This increase in G*, 
although small, nevertheless requires an exponentially long time (in the initial value of -j^f) 
due to the exponentially small nucleation rate. In this exponentially long time, the largest 
clusters grow to an exponentially large size, so large that their growth is diffusion limited. 

In diffusion limited growth, the number n of monomers in a cluster grows at a rate 
proportional to linear size, so h is proportional to ni From the above physical framework, 
creation era scaling units [t] and [n] of time and cluster-size are derived: 

[t}oce^^ F J [njoce^ 7 ^. (1.1) 
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In particular [n] oc This is connected to diffusion limited growth: The ODE h oc ns 

3 

has a solution proportional to t~s. It is perhaps natural that [t] and [n] are exponentially 
large in -j^f, but what is distinctive of the physical model are the pre-factors | and | in 
the exponents in (11. II) . They have in fact been seen before: In a short note, Marchenko 4| 
derived the same characteristic time and cluster-size, with the signature | and |. 

Once physical scaling units are established, an asymptotic theory for the evaluation of 
the cluster size distribution during the creation era readily follows. The mathematical form 
of this theory is the creation signaling problem briefly described in the abstract. Its solution 
is straightforward. The main result is a determination of the nucleation rate as a function 
of time in scale-free form. 

This paper is describes the beginning of the aggregation process as a whole. The long- 
time limit of the creation era provides effective initial conditions for successive stages of the 
aggregation process, eventually making contact with the classic theory of coarsening (aka 
ripening) due to Lifshitz and Slyozov. This subsequent work has been carried out, and is 
being prepared for publication. 



II. THE PHYSICAL MODEL OF NUCLEATION AND GROWTH 

The physical model used here is based on several assumptions: 

Al. The initial super-saturation, and monomer chemical-potential are small (and positive). 

A2. Clusters nucleate at the Zeldovich rate, which adjusts to new values of super-saturation 
immediately. 

A3. The total monomer density (including monomers that form clusters) is conserved. 
A4. The cluster growth is diffusion limited. 
A5. Initially, there are no clusters. 
A6. The temperature, T, is fixed. 

While the growth of clusters is not a deterministic process, the fluctuations are small in the 
size- and time-scales that we investigate here; They can be safely ignored throughout the 
creation era. The evolution of the cluster-size distribution is modeled as an advection PDE, 
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with no diffusion term, in the space of time t, and cluster-size n. The advection velocity 
is the cluster growth rate, and the nucleation rate serves as the boundary condition at size 

71 = 0. 



A. Super-saturation and chemical potential 

Different aspects of the physical model are related via their dependence on the chemical- 
potential and super-saturation. Let f\ denote the density of monomers in the monomer- 
bath, and f s denote the saturation density, that is, the density of monomer that would be in 
equilibrium with an infinite cluster. The super-saturation (p and monomer chemical-potential 
(in units of ksT), r] are defined as 

V=^ L , = log f (2.1) 

Is Js 

Thus, 7] and (p are related by: ip = e v — 1. For 77 <C 1, we have 77 ~ (p. In this section 
we introduce the model using 7/ or <p as appropriate, but in the rest of the paper we use 7/ 
exclusively. 

B. Nucleation rate 

According to Zeldovich jl^ . super-critical clusters are nucleated at a rate per unit volume 
j, given by 

J =tofs X -^-e~^. (2.2) 

Here, f s is the saturation density of the monomer, u is the evaporation rate constant so that 
tun 2 / 3 is the rate at which monomers on the surface of an n— cluster leave it, and a is the 
surface-tension constant, so that fceTcro 2 / 3 is the free-energy associated with the surface of 

3 

an n— cluster. The exponent is the asymptotic approximation of the free-energy barrier, 
k B T- 

C. Growth rate 

As is shown further in this paper, the characteristic size of a cluster during the creation 
era is exponentially large (in G*). This explains our assumption that it is large enough 
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that the clusters' growth is limited by the diffusion of monomers, rather than by surface 
reactions [14]. This is the basis of the Lifshitz-Slyozov formula^ 

n = d(p5 - a ), d= (3(4vr) 2 )^ (Dv* f 8 ) (2.3) 

which describes the growth-rate of an n— cluster. In f |2.3|) . D is the diffusion constant of the 
monomers and v is the volume per particle in the clusters. During the creation era, the first 
term in ( 12. 3ft dominates the second, which is henceforth ignored. However, a balance between 
the two terms reveals the critical size, = (~j , which corresponds to the maximum of 
the free-energy. Ignoring the second term in the formula for n (12.31) . is equivalent to the 
assumption that the characteristic cluster size, [n], is much larger than the critical size n*. 



D. Advection signaling problem 

We approximate the discrete cluster-size densities with a continuous density function 
r(n, t). For small 5n, the density of clusters of size between n and n + 5n is r(n, t) 5n. This 
approximation allows us to write an advection PDE for the cluster-size distribution using the 
growth rate as the advection "velocity". Mathematically, this means that the distribution 
r(n,t) satisfies the PDE 

d t r + dr]d n (n3r) = in n > 0, (2.4) 

where d is defined in ( 12. 3p . The initial conditions we assume are pure monomer, correspond- 
ing to r(n, 0) = for all n > 0. The flux of clusters is evidently drjn^r, and it must tend to 
the nucleation rate as n — > + , giving the effective boundary condition: 

i _ „ 3 

drjn^r — > f2e ^ , as n — > . (2.5) 

Here, VL = fsUj^J^, the prefactor in the Zeldovich formula (12. 2p . The superscript + indicates 
that the limit is taken from above. Readers who are concerned about our use of instead 
of here or in the next subsection, are referred to Appendix O for a brief discussion. 

The advection PDE (12. 4p subject to the effective boundary conditions (12. 5p and the initial 
condition r(n, 0) = constitute the creation signaling problem mentioned in the abstract. 
It determines the evolution of r(n, t) as a functional of rj = r](t). The model is "closed" by 
a determination of rj(t) as a functional of r(n, t) using the conservation assumption. 
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E. Conservation of total monomer density 



The conservation of monomers is expressed (approximately) by 

oo 

f = fi + J nr(n,t)dn. (2.6) 
o 

Here, the total monomer density /, a constant, is the sum of monomer density, fi in the 
bath, and J °° nr(n, t) dn, which approximates the density of monomers in clusters. Inserting 
the relation fx ~ (1 + r))f s , which follows from (12 .ip in the limit 77 <C 1, into (12.61) we find 

00 

f = (l +v )f s + J nr{n,t)dn. (2.7) 
'0 

For a full derivation and further discussions on these models, we refer the readers to Wu's 
review article [3] and references therein. 

III. ASYMPTOTIC SOLUTION OF THE CREATION SIGNALING PROBLEM 

The equations describing the creation era advection signaling problem are gathered to- 
gether for quick reference: 

d t r + d n (dnn^r) = in n > 0, (3.1) 
r(n,0) = 0, for n > 0, (3.2) 

1 _ CT 3 

dr/n'^r — * j = Qe ^ as n — > , (3-3) 

00 

r)(t) = T)(0) - \- I nr(n, t) dn. (3.4) 





We assume that the initial super-saturation is small, and take e = 7](0) C 1 as the small 
parameter of the asymptotics. 

A. Dominant balance scalings 

In the limit < r\ <C 1, the nucleation rate in (13.31) is highly sensitive to small changes in 
the chemical potential 77. Hence, we work with the change in chemical potential, 

Srj = r] — r/(0) = i] — e, 



and expect the super-saturation n to remain close to its original value 77(0) throughout the 
creation era. We find scaling units [t] , [n] , [r], and [Srj\ of the variables t, n, r, and 5rj from 
dominant balances applied to equations fl3.1H3.4l) . 

The dominant balance associated with the advection PDE (13.1 p is 

frv\ ick[n]i) - (3 - 5) 

Integrating (13. ip from n = to 00 and using boundary condition (13.31) we find 

CO _ a 3 

4 J r dn = fie , with corresponding dominant balance 


L J L J -. Qe~^. (3.6) 



00 

The conservation equation ( 13. 4ft can be written as Si] = ^ J nr(n, t) dn, with corresponding 

/s 

dominant balance 

N = |[r][n] 2 . (3.7) 

In addition to the three "straightforward" dominant balances (l3.5H3.7p which follow directly 
from ( I3.1H3.4I) . there is one that quantifies the change Sn in chemical potential required to 
"shut off" nucleation. For 5rj <C rj, the relative change in nucleation rate j in ( 13.31) that 
results from a change 5n in rj is U- ~ ^Srj, so the scaling unit [5rj] of £77 that corresponds to 
a significant change in the nucleation rate is 

N = 4 (3-8) 



This is consistent with our expectation that it is small relative to rj. We substitute (13. 8p 



into (13. 7p . and then solve (I3.5H3.7I) for the scaling units [t] , [n] , and [r] of time, cluster size, 
and cluster-size density: 



i/7. 



2 

5 2 G, 
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Recalling the definitions of d and f2, these become 



[t] = (8n)-r, jet e^ fc sT [D 3 vf^u 2 ) « 



[n] = (tt^2^\/3) 



3 _G* 
g 5 k fl r J 



[r] = (3 • 2 1 V)- 1 / 5 



(3.9) 
(3.10) 

(3.11) 



In the exponents, = ksTj^ is the initial free energy barrier. 



B. The Reduced Kinetics 

We nondimensionalize the creation signaling problem (13.11 - 13. 41) using the units in (13. 8f - 
13. lip . In the limit e — > 0, the reduced equations are: 

d t r + d n (nh) = 0, in n > 0, (3.12) 
n^r^e' 5 ' 7 , as n -> + , (3.13) 

oo 

Sr] = — nrdn. (3-14) 



o 

The initial condition is 

r(n, 0)=0. (3.15) 

The reduced signaling problem (I3.12H3TT5I) is transformed into an integral equation for 8r], 
which is solved numerically. The density r(n,t) is subsequently recovered from 5rj(t). 
The flux of super-critical clusters, n^r, is constant along the characteristics 

^ - T )) a '*) : ^ r }' ( 3 - ie ) 

of the PDE (I3.12p . The characteristics C T can be seen in Fig. [TJ in which the (horizontal) 
density of the characteristics at each point (n, t) is proportional to the density of clusters of 
size n at time t. The curves in ( 13 .161) describe the world-lines of clusters that nucleate at 
time t = t. The region below the thick line in Fig. [1] corresponds to t < |nt . In this region, 
there are no cluster world lines and r(n, t) = 0. 




2 4 6 8 10 



FIG. 1: The characteristics of the reduced PDE. The flux mr is constant along the characteristics. 
The density of the curves is proportional to the solution r(n, t). 



For a known 8r)(t), the solution r(n, t) that has a constant flux n\r along characteristics, 
and satisfies the BC (13.131) . is 

r(n, t) = { % (3.17) 



0. 



3 2 
£ > -77,3, 

~ 2 ' 

3 2 

< £ < -«3. 

2 



An integral equation for 5rj(t) is found by substituting (13.171) for r(n, t) in the conservation 
identity (I3T41) : 

t 3 

r /n \ 5 

(3.18) 




In (I3.18P the variable of integration has been changed from n to r = £ — |n 2//3 . 



C. Creation Transition and Physical Predictions 

We solve (I3.18P numerically. A short discussion of the method and numerical result can 
be found in Appendix |A] The nucleation rate 

j = e *i® i s 

calculated and plotted in Fig. [21 
At time t = 2, j is about a third its original value, and at t = 5 it has effectively vanished. 
The distribution r(n,t) of cluster sizes is recovered from 8rj(t) via (13.171) . Figure [3] shows r 
vs. n for an increasing sequence of t. 

From the numerical solution for j we find the total density of clusters generated during 
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i(t) 




FIG. 2: The Zeldovich nucleation rate j as a function of time. After an 0(1) time, the super- 
saturation decreases slightly, and the nucleation rate gets turned off. 




FIG. 3: The density of cluster sizes, r(n, t), for various values of t. 



the creation era. This density is given by the integral J rdn. Using the same change of 

o 

variables from n to r as in (13.171) and PDE (13.121) . we convert the integral of r into an 
integral of j = e Sv , 

oo oo 

R = Jrdn = J j(r)dr. (3.19) 
o o 

The value of R, based on the numerical approximation to j(t), is 



1.7109. 



Converting back to original physical units, the total density of clusters produced during the 
creation era is 



R\r]\n] 



3^R 



2i0 7Tio {f s Dv3(j 



I 3 
2 



(3.20) 
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D. The Long-Term Behavior 



We discuss the long-term behavior of the cluster-size distribution, r(n, t), is determined 
by the long-term behavior of j through Eq. (13.171) . For t large relative to the "main support" 

of j = e Sr) , we approximate the integral equation (13. 18j) by 

3 i r°° 

-5ri~ (ft) 3 fl+(ft) 3 / re^ (T) dr. (3.21) 

Jo 

Fig. [2] shows the graph of j(t) together with this approximation. The second moment of 
j(r) is found to be 

oo 

Tj{T)&Ttt 1.7773. (3.22) 

Since j(t) diminishes to a small fraction of its original value after t — 5, the main support 
of r(-, t) is located in the range of n where < t — < 5. This can be seen most clearly 
from (13.171) . The upper end of this interval is sharp: At time t there are no cluster sizes 
greater than (§t)^, but the lower end is "soft" since although j(t) is small at t — 5, it is not 
exactly zero for any t > 0, however the super-exponential decay of j(t) allows us to treat the 
support of j (and therefore of r) as being well-defined despite this softness. By comparing 
the width of the distribution to the size of the largest cluster N, 

N(t) = fct\ 2 , (3.23) 

we show that outside of a (relatively) narrow region the distribution is exponentially small. 
In this sense we can say that the support of r is concentrated in a narrow interval, and the 
distribution is asymptotically monodisperse. Describing this narrow region is our final task. 
We define the distance between n and the size of the largest cluster iV: 

Sn = N - n. (3.24) 

The support of r is concentrated at values of 5n for which < t — |(iV — 5n)i < 5. Thus, 
to first order in the support of r is concentrated where < < 5. In other words, the 
width of the distribution grows like t? but the size of the largest cluster (i.e., the location 

3 

of the distribution) grows like t?, so the relative width of the distribution is shrinking. The 
long-term description of the cluster-size density is given by 

r(n , t) = { N --> M ■ ° < ^ < 5 (3 , 5) 

0, otherwise. 
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Here, N and 5n are defined in (I3.23j) and (13.241) . The constant 5 is used to describe the 
(soft) upper bound of the support of j(t). 



IV. CONCLUSIONS 

The exponential dependences of the characteristic time and cluster size upon the initial 
free energy barrier are based on order of magnitude balances applied to simple, essentially 
classical kinetics. They should be the most robust results of this paper. 

We mention a biophysics application: Hydrophobic proteins have been "crystallized" into 
periodic arrays for x-ray analysis. They are implanted into a space filling cubic phase of 
lipid bilayer, which acts as the "solvent" . Nucleation of "protein crystals" is observed. The 
size of the cubic phase unit cell was manipulated, and it was observed that the characteristic 
size of the crystals, and their time of formation decrease as the size of the unit cell increases. 



lie 
21 



Using the proposed exponential dependences of crystal size and time of formation upon t 
free energy barrier, the second author, together with M. Grabe, G. Oster, and P. Nollert 
quantified the decrease of free energy barrier with increasing size of unit cell. The results 
are consistent with an independent estimate of the free energy barrier based upon the elastic 
energy of embedding proteins in the bilayer. 

Certain detailed results of the current paper are expected to be less robust, and should be 
regarded as documenting the predictions of (oversimplified) classical kinetics. In particular, 
the emerging distribution of cluster sizes after nucleation becomes conspicuously narrow, so 
the sizes of clusters are much closer to one uniform size than is observed in experiments. One 



proposal [TjJ is that a broader distribution results if in fact we observe the superposition of 
distributions in a spatially inhomogeneous nucleation process. 

The sharp front of the cluster size distribution at the largest cluster size as predicted by 
our model is also expected to be a casualty of any comparison with reality. Many effects could 
regularize it. In particular, B. Niethammer and J.J.L. Velasquez [6| formulated a diffusion- 
like correction to the advection PDE of the LS model, based on screening fluctuations in the 
local super-saturation seen about clusters. 

The extreme sensitivity of the nucleation rate upon small changes in the chemical po- 
tential rj of monomers was exploited in our asymptotic solution of the creation transient. 
But this extreme sensitivity is a potential source of difficulties as well. For instance, most 
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works to date, such as Penrose 



7j, Niethammer and Velasquez use the approximation 
to the conservation of particles, which says: "Density of monomers plus density of 'large' 
clusters equals total particle density." In our work, this simple approximation to particle 
conservation is retained, so as to not distract from the main results. However, even small 
corrections to the conservation law and a consequential small change in rj can be amplified 
to large (relative) corrections to the nucleation rate, which is exponentially small as 77 goes 



to zero. Specifically, Wu in 



12| mentions the quasi-equilibrium distribution of sub-critical 



clusters ('embryos' in his terminology) with 1 < n < n*. We propose that the inclusion of 
the embryos in the conservation of particles is one of those small corrections that lead to 
significant changes in predictions of the nucleation rate. 
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Appendix 

APPENDIX A: NUMERICAL SOLUTION OF INTEGRAL EQUATION 

The analysis of the nucleation era requires the solution of the integral equation, 

Sv(t) = ~ I (| (t - t)) 2 e***> dr. (A.l) 


We find anapproximate solution of this equation on a set of equally spaced times t n on 
the interval [0, 15]. The value of Sr](t n+ i) is calculated from the integral of 5r] up to time 
t n+ \ using the trapezoidal rule approximation. This is not implicit as the t — r term in the 
integrand of (lA.lj) vanishes at r = t. The order of accuracy of the method is found to be 
I (see Fig. |4|). This fractional order is probably due to the cusp in the integrand at t — r, 
which adds an error of 0(At) 3 / 2 to the integral. The convergence of the error can be seen 
in Fig. H] where the (natural) log of the errors (of both j and R) are plotted against the log 
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-6 

o 

aj _8 




6 8 10 12 

log number of points 



FIG. 4: The errors in j and R, together with a line of slope —3/2 in a log-log plot. The x-axis is 
the log (base 2) of the number of points used, and the y-axis is the (natural) log of the error. 

(base 2) of the number of points in the segment [0, 15]. The most accurate result we have 
obtained — with 2 13 points — for the resulting density of clusters is R = 1.7109162 ± 3. 

APPENDIX B: VALIDITY OF THE QUASI-STEADY-STATE ASSUMPTION 



One of the assumptions used in this paper is that the nucleation rate corresponds to the 
Zeldovich rate with the instantaneous super-saturation. This "quasi-steady-state" (QSS) 
assumption has been previously studied by Shneidman [9|. It was shown that the QSS 
assumption is valid if 



(B.l) 



where 7 is given by 



7 



d 

-T — 



dt \ ksT 



dh 
dn 



(B.2) 



In the case of the creation era, we have shown that the characteristic time [t] is exponentially 
long, as given by (13.91) . Therefore the derivative of the free-energy barrier is exponentially 
small. In flB.2j) . r is a parameter that depends on the physical model for the growth of 
small clusters. For example, the Becker-Doring model, has r = Thus, we see that the 
exponent 7 is itself exponentially close to zero (due to the time derivative), and the LHS in 
fIB.ll) is extremely close to 1. This validates, a-posteriori, the QSS assumption. 
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APPENDIX C: CRITICAL SIZE n* IS NOT EQUAL TO ZERO 



The effective boundary condition (12.51) results from a heuristic asymptotic matching in 
the small super-saturation limit rj — > 0: The Zeldovich rate on the RHS is a quasi-static 
discrete flux which measures the net rate of creation of n + 1— clusters from n— clusters per 

3 

unit volume, valid for n on the order of critical size n* = The LHS is a continuum 
approximation to the flux based on diffusion limited growth. Presumably, diffusion limited 
growth is valid for some range of "large" cluster sizes with n ^> n*. Of course we assume that 

3 T 3 

the characteristic size [n] <x e 5 <C is in this range. Moreover, we assume that diffusion 
limited growth holds asymptotically for clusters of an "intermediate" size n + = n + (rj) so 
<C n + <C [n]. Consider clusters with size n in the interval n* < n < n + . Under 
quasi-static conditions, we expect that influx through the n = end at the Zeldovich rate 
balances the "diffusive limited" outflux through the n = n + end, so that 

_ - 3 i 

fie = dr]n\r(n + ,t). (C.l) 

Since n + /[n] ^ as rj — > 0, we obtain (12.51) as the formal limit of (1C.1I) as r] — > 0. 

The appearance of zero as the lower limit of the integral J °° nr(n, t) dn in (2.6) is also 
part of the small super-saturation limit rj — > 0. The exact density of monomers in clusters 
is the discrete sum Y^=2 n fni where f n is the density of n— clusters. The standard idea 
behind replacing this sum with an integral is that the f n are the values of a smooth function 
at integer arguments, whose characteristic scale of the independent variable is much larger 
than unity. That is almost what we have. We approximate f n by r(n, t) whose characteristic 

3 CT » 

scale of n is [n] oc e 5 ^? — > oo as rj — > 0. Of course the continuum approximation to f n 
breaks down for some range of n with n <C [n], so in addition we are assuming that the 
contributions to the sum Yl^=o n f n fr° m this n [n] range are negligible as rj — > 0. 
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